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^ ■ Abstract 



Suzaku deep observations have discovered two highly significant nonthermal X-ray 
sources, Suzaku J1804-2142 (Src 1) and Suzaku J1804-2140 (Src 2), positionally 
coincident with the unidentified TeV 7-ray source, HESS J1804— 216. The X-ray 
sources are not time variable and show no counterpart in other wavebands, except 
for the TeV source. Src 1 is unresolved at Suzaku spatial resolution, whereas Src 2 
is extended or composed of multiple sources. The X-ray spectra are highly absorbed, 
hard, and featureless, and are well fitted by absorbed power-law models with best- 
fit photon indices and absorption columns of — 0.3loJ and 0.21^2 x 10 22 cm -2 for 
Src 1, and 1.7±H and x 10 23 cm for Src 2. The measured X-ray absorption 

to the latter source is significantly larger than the total Galactic neutral hydrogen 
column in that direction. The unabsorbed 2-10 keV band luminosities are 7.5 x 
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10 32 (g?/5 kpc) 2 ergs s -1 (Src 1) and 1.3 x 10 33 (g?/5 kpc) 2 ergs s" 1 (Src 2), where d is the 
source distance. Among the handful of TeV sources with known X-ray counterparts, 
HESS J1804— 216 has the largest ratio of TeV 7-ray to hard X-ray fluxes. We discuss 
the nature of the emission and propose the Suzaku sources as plausible counterparts 
to the TeV source, although further observations are necessary to confirm this. 

Key words: acceleration of particles — ISM: individual (HESS J1804— 216) - 
X-rays: ISM — X-rays: individual (Suzaku J1804-2142, Suzaku J1804-2140) 

1. Introduction 

The origin and acceleration mechanism of cosmic rays has remained a long-standing 
problem since their discovery nearly 100 years ago. The most widely accepted acceleration 
mechanism, at least for cosmic rays up to 10 15 ' 5 eV (the knee energy), is diffusive shock ac- 
celeration at the rims of young supernova remnants (SNRs) (e.g., Bell 1978; Blandford & 
Ostriker 1978). This hypothesis was greatly strengthened by the discovery of localized regions 
of synchrotron X-ray emission from the rims of SNR SN 1006 (Koyama et al. 1995), which 
indicated that SNRs were able to accelerate electrons up to ~TeV energies. Since then, several 
other SNRs have been identified as electron accelerators (e.g., Koyama et al. 1997; Slane et 
al. 2001; Bamba et al. 2000; Bamba et al. 2001) and numerous studies have indicated that the 
acceleration process is very efficient (e.g., Vink & Laming 2003; Bamba et al. 2003; Bamba et 
al. 2005; Volk et al. 2005). Many papers regard RX J1713— 3946 as a proto-type of cosmic-ray 
accelerating SNRs, since it is bright in both X-ray and TeV bands (Uchiyama et al. 2003). To 
date the best evidence for the acceleration of hadronic cosmic rays comes from observations of 
the Tycho SNR (Warren et al. 2005). 

Recently, a TeV 7-ray telescope, High Energy Stereoscopic System (HESS), has discov- 
ered some new sources in the Galactic plane (Aharonian et al. 2005a; Aharonian et al. 2006). 
Some of these sources are extended on scales of ~ 0.°1, suggesting Galactic counterparts like 
pulsar wind nebulae (PWNe), SNRs, star- forming regions, and so on. However, their nature 
is still unclear: mainly because many of the new sources have no counterparts in other wave- 
bands. They have been sometimes referred to as "dark particle accelerators" . Deep follow-up 
observations are crucial to elucidate their nature and hard X-ray observations are a valuable 
tool for this, since above a few keV we do not have to contend with absorption by the Galactic 
plane. The X-ray satellite Suzaku has large effective area and low background (Mitsuda et 
al. 2006), making it an excellent choice for such follow-up studies. 

HESS J1804— 216 is one of the brightest of the unidentified (unID) TeV sources 
(Aharonian et al. 2005a). Its size (0.°2) and photon index (2.72 ±0.06) makes it one of the 
largest and softest sources. Aharonian et al. (2006) proposed two possible counterparts: the ra- 
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dio SNR G8.7-0.1 (W30) and a young Vela-like pulsar PSR J1803-2137, although neither is a 
definitive identification. The HESS source overlaps only the southwestern portion of G8.7— 0.1, 
which is about 3/4° in diameter. The distance is estimated to be ~5 kpc based on associating 
the SNR with nearby H II regions (Kassim & Weiler 1990; Finley & Oegelman 1994). From CO 
observations, the surrounding region is known to be associated with molecular gas where mas- 
sive star formation takes place (Ojeda-May et al. 2002). From the observed dispersion measure, 
PSR J1803— 2137 is located at a distance of ~5.3 kpc (Kassim & Weiler 1990). Both candidates 
were detected by ROSAT (Finley & Oegelman 1994); the soft X-ray emission of G8.7— 0.1 is 
only apparent along the northeastern radio rim. Very recently, Brogan et al. (2006) reported 
the discovery of a new SNR (G8.31— 0.9) in this region. Swift independently discovered some 
X-ray sources in this region recently (Landi et al. 2006), but the statistics was quite limited 
due to the short observation. It is interesting to note that the AGASA and SUGAR groups 
have reported an excess of EeV cosmic rays from the general direction (i.e., within ~20°) to- 
ward HESS J1804-216 (Hayashida et al. 1999; Bellido et al. 2001), although this result has 
not been confirmed by more sensitive observations by the Pierre Auger Observatory (Letessier- 
Selvon al. 2005). The evidence that HESS J1804— 216 is a proton accelerator therefore remains 
ambiguous (Fatuzzo et al. 2006). 

In this paper, we report the first results of Suzaku follow-up observations of 
HESS 1804—216. §2 introduces the observations, our analysis is explained in §3 and we discuss 
our findings in §4. 

2. Observations and Data Reduction 

We observed a source and background region with Suzaku (Mitsuda et al. 2006), which 
has two sets of active instruments: four X-ray Imaging Spectrometers (XIS, Koyama et al. 2006) 
each at the focus of an X-Ray Telescope (XRTs, Serlemitsos et al. 2006) and a separate Hard X- 
ray Detector (Takahashi et al. 2006; Kokubun et al. 2006). Here, we present the analysis of XIS 
data only. For this study we observed both a source region (centered on the published position 
of HESS J1804— 216 and covering the brightest part of the TeV source) and a background region 
some 30' to the southwest that avoided the TeV source (see Figure 1). Observational details 
are given in Table 1. 

The XIS was operated in normal full-frame clocking mode in both pointings. The data 
reduction and analysis were made using HE AD AS software version 6.0.6. We used version 
0.7 processed Suzaku data, which is an internal processing applied to the Suzaku data during 
the science working group phase for the purpose of establishing the detector calibration as 
quickly as possible (for the detail, see e.g., Fujimoto et al. 2006). The XIS pulse height data for 
each X-ray event were converted to pulse invariant channels using the xispi software version 
2005-12-26 and charge transfer inefficiency parameters from 2006-05-24. We filtered out data 
obtained during passages through the South Atlantic Anomaly or with elevation angle to the 
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Earth's limb below 5°. Events with ASCA grades of 0, 2, 3, 4, and 6 were retained. The total 
available exposures are ~40 ks and ~43 ks for source and background, respectively. 

3. Results 

3.1. Images 

Figure 2 shows the XIS images in the 0.5-2.0 keV (left) and 2.0-7.0 keV (right) bands. All 
four XIS images were combined, corrected for vignetting, and smoothed with a ~20" (sigma) 
Gaussian. In the soft X-ray band, we see no significant source at the best-fit position of 
HESS J1804— 216 (large cross). On the other hand, there are two discrete sources in the 
hard band image (the right panel): at positions (epoch J2000) of (18 h 04 m 44 s , -21 d 42 m 42 s ) and 
(18 h 04 m 34 s ,-21 d 40 m 20 s ). We designate them as Suzaku J1804-2142 and Suzaku J1804-2140. 
Note that the systematic position uncertainty is about 1 arcmin. Hereafter, we will refer to 
these as Src 1 and Src 2, respectively. In order to examine their significance, we extracted 
source photons from a circle of 1.5 arcmin radius for Src 1 and 1.8 arcmin radius for Src 2 
(see circles in right panel of Figure 2). The total 2.0-7.0 keV count rate was 3.7xl0 -2 cnt s _1 
and 5.3xl0~ 2 cnt s" 1 , respectively. The background counts in the same energy range was 
estimated from the background observation as 2.8 xlO -3 cnts s~ 1 arcmin -2 . The significance 
of Src 1 and Src 2 are, then, 12.4<r and 14.9<r, respectively. The Galactic ridge emission is 
known to have spatial fluctuations of about 30% in intensity (Yamauchi et al. 1996; Kaneda 
et al. 1997). Assuming a background rate that is 30% higher than our nominal one, does not 
change the source detection significances critically. Data obtained at low geomagnetic cut-off 
rigidity (COR) values tends to have higher background, so different ranges of COR values in the 
source and background pointings could introduce a difference in background level. We checked 
the importance of this by heavily filtering on COR (allowing only values greater than 10 GV) 
and found that we could still detect the sources significantly even after such a drastic reduction 
of the data. 

The spatial extent of the sources were also examined. The half power diameter of the 
XRT is about 2 arcmin and is almost independent of X-ray energy (Serlemitsos et al. 2006). 
We compared the counts within the circle of 1 arcmin radius with the counts from within a 
radius of 1.5 arcmin (Src 1) and 1.8 arcmin (Src 2). According to XRT calibration data the 
ratio of counts from a point source for these radii should be 71% and 63%. The background 
estimation was done in the same way as discussed above. We found that about 75% (Src 1) 
and 42% (Src 2) of the counts are concentrated within the 1 arcmin-radius circle. Therefore, 
we conclude that Src 1 is point-like or compact compared to the spatial resolution of Suzaku, 
whereas Src 2 is significantly extended or multiple. 

Src 2 is spatially closest to the published position of HESS J1804— 216 (see Figure 2), 
while both Src 1 and Src 2 are within the brightest parts of the TeV source. No counterparts 
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were found in the SIMBAD database, catalogs of H II regions (Ojeda-May et al. 2002; Kassim 
& Weiler 1990) and ultra-compact H II regions (Bronfman et al. 1996). Figure 1 shows the 
TeV 7-ray image of HESS J1804— 216 (grayscale) overlaid with contours of the XIS 2.0-7.0 keV 
(white) and soft X-ray (ROSAT) emission (black). Both ROSAT and ASCA imaged the Suzaku 
XIS field of view, but failed to detect emission from either Src 1 or 2, due to low detection 
efficiency (ROSAT) in the hard X-ray band and strong stray light from a nearby bright source 
(ASCA). the Suzaku sources are quite compact compared to the HESS source. Note that 
PSR 1803—2137 and the previously detected soft X-ray emission from the northeastern shell of 
G8. 7-0.1 fell outside the XIS field of view (Figure 1). 

We extracted light curves for the new X-ray sources, rejecting data obtained in low 
COR regions (<10 GV), and applied a Kolmogorov-Smirnov test for time variability. No 
significant variability was detected for both sources during the observation, with the probability 
of constancy of 0.38 (Src 1) and 0.49 (Src 2). 

3. 2. Spectra 

Source spectra were extracted from within circular regions of radius of 1.5 arcmin (Src 1) 
and 1.8 arcmin (Src 2). Background spectra were obtained from the background observation 
in order to account for the local Galactic ridge emission. The background spectrum was accu- 
mulated from within a 3 arcmin radius circle centered on the aim-point of this second pointing. 
The background rate in different observations can differ due to differences in the average COR. 
However our observations were done within one day of each other and the attitude of the satel- 
lite was very nearly the same. The average COR differs by only a few percent between the 
source and background observations, so we have ignored any effect of the COR range on the 
background. 

As shown in Figure 3, the background-subtracted spectra are highly absorbed and fea- 
tureless, so we fitted them with absorbed power-law spectral models. The absorption model 
used the cross sections of Morrison & McCammon (1983) with solar abundance (Anders & 
Grevesse 1989). We made auxiliary files for the effective area using xissimarf gen (Ishisaki et 
al. 2006), which includes the effect of contamination (Koyama et al. 2006). Figure 3 and Table 2 
shows our best-fit models and parameters within 90% errors. In each case the best-fit result 
was acceptable statistically with reduced x 2 values of 47.1/46 for Src 1 and 34.6/38 for Src 2. 
It was also possible to obtain acceptable fits using a non-equilibrium ionization thermal plasma 
model (Borkowski et al. 2001); however the derived temperatures, 80 (>21) keV for Src 1 and 
14 (>2) keV for Src 2, were unrealistically high. We conclude that a thermal interpretation for 
these spectra is unlikely and favor a nonthermal interpretation instead. 

Given the large mis-match in sizes between the new X-ray sources and HESS J1804— 216 
we also investigated how much diffuse hard X-ray emission there could be extending over the 
entire XIS field of view. We extracted photons from the entire field of the source pointing, 
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excluding Src 1 and Src 2. A background spectrum was made using all the photons in the 
background pointing. There was no significant excess in the source pointing compared to the 
background over the 2.0-10.0 keV band. Assuming a power-law model with fixed photon index 
(T = 2) and high absorption (Nh = 10 23 cm" 2 ), we obtain an absorbed flux upper limit (within 
a 90% error) of 5 x 10 -13 ergs s _1 cm -2 . 

4. Discussion 

The new hard X-ray sources we have successfully detected with Suzaku are now the near- 
est potential counterparts to HESS J1804— 216. Although the morphology of TeV 7-ray emis- 
sion and hard X-rays are quite different, this is also the situation for the radio SNR G8.7— 0.1 
and the soft X-ray emission. None of the possible counterparts is a morphological match to the 
HESS source. Nevertheless the new X-ray sources provide interesting new information on the 
nature of HESS J1804— 216 and the environment. An independent report (Landi et al. 2006) 
notices new X-ray sources in this region using Swift. The position and the flux of the source 
nearest to HESS J1804— 216 is consistent with those of Src 2, indicating that we successfully 
confirmed their detection. We estimated the chance coincidence probability of unrelated X-ray 
sources for HESS J1804— 216. With the logiV-logS' relation of X-ray sources in the 2.0-10.0 keV 
band (Sugizaki et al. 2001; Yamauchi et al. 2002; Ebisawa et al. 2005), we found that the ex- 
pected number of X-ray sources is only 4-9 xlO -3 in the elliptical region with radii of 0.016 and 
0.018 degrees (the position uncertainty of the HESS source; Aharonian et al. 2006). Therefore, 
we concluded that they are physically associated with HESS J1804— 216. In the following we 
discuss Src 1 and Src 2 separately. 

The best-fit absorbing column for Src 1, 0.21^2 x 10 22 cm -2 , is consistent with the 
Galactic neutral hydrogen column in this direction of 1.5 x 10 22 cm -2 (Dickey & Lockman 1990). 
The unabsorbed luminosity is 7.5 x 10 32 (<i/5 kpc) 2 ergs s _1 in the 2.0-10.0 keV band, where 
d is the source distance. We quote results for a nominal value of 5 kpc although in fact we 
have no independent constraints on the distances to the new X-ray sources. Src 1 has a very 
flat spectrum, flatter than typical for young pulsars and their associated wind nebulae (PWN) 
(typical T values of 1.5 to 2.5; Possenti et al. 2002) or shell-like SNRs emitting synchrotron 
X-rays (typical V values of 2.3 to 3). Uchiyama et al. (2002) have reported a region near SNR 
G78. 2+2.1 (Gamma-Cygni) with a photon index of 0.8-1.5, that might be comparable to Src 1. 
High mass X-ray binaries also have hard photon indices (r of 0.5-2.5; e.g. Muno et al. 2003). 
So Src 1 might be this type of source. 

In contrast Src 2 has a best-fit absorbing column that is about an order-of-magnitude 
higher than the expected Galactic column. This is consistent with Src 2 being embed- 
ded in dense gas, although there are no cataloged molecular clouds in this direction. The 
X-ray and radio counterpart to HESS J1813— 178 (Brogan et al. 2005) also appears to 
show evidence for excess absorption. The unabsorbed 2-10 keV band luminosity of Src 2, 
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1.3 x 10 33 (<i/5 kpc) 2 ergs s _1 , best fit photon index (see Table 2), and extended nature are all 
consistent with either SNR origin: a PWN or shell-like synchrotron emission. 

Their origin is still unknown, but there is some interesting scenario for their nature; old 
SNRs colliding with molecular clouds. Yamazaki et al. (2006) argue that an old SNR shock 
(age of 10 5 yr) interacting with a giant molecular cloud can emit hard nonthermal X-rays 
and strong TeV 7-rays from shock accelerated protons through ir° decay (see also Fatuzzo et 
al. 2006). They can explain large TeV 7-ray to X-ray flux ratios (factors of ~100). Src 2 might 
represent the X-ray component to the HESS source under this scenario, with the large excess 
absorption as a key piece of evidence. 

Comparing images in other wavelengths also suggests the SNR origin. Soft X-ray emis- 
sion and HESS source are both encompassed by the radio contours of G8.7— 0.1, suggesting 
a common origin, although they do not overlap, On the other hand, the heavily absorbed 
Suzaku-discovered hard X-ray sources are located right at the position of peak HESS emission. 
The anti-correlation of HESS and ROSAT emission is likely a result of a dense region near the 
SNR which absorbs the soft X-rays, while at the same time enhancing the production of TeV 
emission from ir° decay. The larger size of the TeV source compared to the hard X-ray source 
might be due to the much larger diffusion length of protons (which provide the TeV 7-rays) 
compared to primary electrons (which emit the hard X-rays). 

HESS J1804— 216 has the large ratio of TeV 7-ray to hard X-ray fluxes, which should 
be an important key to solve its nature. Table 3 summarizes X-ray and TeV 7-rays fluxes 
for a sample of TeV 7-ray emitters. Generally speaking, PWNe have rather small values of 
Frev/Fx, whereas SNRs tend to be brighter in TeV 7-rays. For HESS J1804— 216, we assumed 
three counterparts; Src 1, 2, and the upper limit of the diffuse emission. We extracted the diffuse 
emission from the entire region of XIS FOV, whereas the TeV emission comes from wider region 
(er = 0.2 degree; Aharonian et al. 2006). Therefore, we re-normalized the flux by assuming that 
X-ray comes uniformly from the same region of TeV 7-rays. The flux in the TeV band is 
the integrated one from the entire region of the HESS source. The ratio of HESS J1804— 216 
is significantly larger than that of PWNe (Willingale et al. 2001; Aharonian et al. 2004a), 
and known shell-like SNRs emitting synchrotron X-rays (Ozaki & Koyama 1998; Aharonian et 
al. 2005d, for example). In fact, the F^ev/Fx value is the largest among the identified objects, 
even if we include the upper-limit to the total hard X-ray diffuse emission from the entire XIS 
field of view The value is consistent with that of HESS J1303— 631, one of the new class of 
"dark particle accelerators" (Mukherjee & Halpern 2005; Aharonian et al. 2005f, see Table 3). 

Under the assumption that the TeV 7-rays are emitted via inverse Compton on relativis- 
tic electrons, higher values of F^ev I Fx suggest weaker magnetic fields, all other things being 
equal. Therefore, the large FxeV / 'Fx of HESS J1804— 216 suggest that the magnetic field in 
this source is weaker than typical in traditional PWNe and SNRs, or inverse Compton origin 
is unlikely. 
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5. Summary 



Thanks to the large effective area and low background of the XIS onboard Suzaku, 
we have discovered two discrete X-ray sources positionally coincide with a HESS unidentified 
source, HESS J1804— 216. The new X-ray sources have no counterpart in other wavebands, 
except for the TeV source. Their X-ray spectra are hard and deeply absorbed. The absorption 
of Src 2 significantly exceeds the Galactic value implying that it is deeply embedded in dense gas, 
which gives strong support to the idea that the TeV 7-rays arise from tt° decay. Additionally 
the Fxev/Fx value is the largest among known sources. HESS J1804— 216 an important object 
for trying to untangle the nature of a "dark particle accelerator." 

Much remains to be done to confirm the nature of these new sources, especially in the 
radio band, which should be able to address the SNR origin. X-ray and TeV 7-ray observations 
with better spatial resolution are also needed and high time resolution will help to judge whether 
there is a pulsar. 
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Table 1. Observation Log. 




source 


background 


Date (yyyy/mm/dd) 


2006/04/06 


2006/04/07 


Aim point (RA, Dec.) 


(271.1713, -21.6759) 


(270.9600, -22.0244) 


Exposure (ks) 


40 


43 
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Fig. 1. HESS image (gray scale; courtesy of W. Hoffman and S. Funk) overlaid XIS 2.0-7.0 keV im- 
age (white contour), with J2000 coordinates. The black contour and circle represents ROSAT image of 
G8. 7-0.1 (Finley & Oegclman 1994) and PSR 1803-2137. All images arc in the linear scale. The black 
boxes show the field of views of XIS. 
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Fig. 2. XIS images of HESS J1804-216 region in the 0.5-2.0 keV band (left) and 2.0-7.0 keV band 
(right), with J2000 coordinates. The scales are linear. Each picture is smoothed with ~20 arcsec scale. 
The thick crosses in both panels represents the best-fit position of HESS J1804— 216. The source regions 
used for the timing and spectral analysis are shown by thin solid lines in the right panel. 
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Table 2. Best-fit parameters of spectral fittings. t 



Src 1 Src 2 



T -0.3 (-0.8-0.2) 1.7 (0.7-3.1) 

./Vh [10 22 cm- 2 ] 0.2 (<2.2) 11 (5.3-21) 

intrinsic F 2 - W kcV* 2.5 (2.1-2.9) 4.3 (3.2-8.3) 

observed -F2-10 kcV* 2.5 2.4 

Parentheses indicate single parameter 90% confidence regions. 
+ : In the unit of 10 _13 ergs s _1 cm -2 . 



Table 3. Flux ratio between X-ray and TeV 7-ray bands. 



F x a F TeV b F TeV /F x Type c Ref. d 

Crab nebula 2.1 x 10 4 56 2.7 x 10~ 3 PWN (1) (2) 

Nebula of PSR 1509 . 32 16 0.5 PWN (3) (4) 

GO.9+0.1 5.8 2.0 0.3 PWN (5) 

SN 1006 NE 19 <2.6 <0.1 SNR (6) (7) 

RX J1713-3946 540 35 0.06 SNR (8) (9) 

HESS J1813-178 7 9 1.3 SNR (10) (11) 

HESS J1303— 631 <6.4 10 >1.6 unID (12) (13) 

HESS J1804-216 . . . . 0.2 (Src 1) 6.5 33 unID this work, (11) 

HESS J1804-216.... 0.4 (Src 2) 6.5 16 unID this work, (11) 

HESS J1804-216 . . . . <0. 8 (diffuse) 6.5 >8 unID this work, (11) 



a : Unabsorbcd flux in the 2-10 keV band in the unit of 10~ 12 ergs s _1 cm~ 2 . 
b : Unabsorbed flux in the 1-10 TeV band in the unit of 10~ 12 ergs s~ J cm~ 2 . 
c : PWN: Pulsar wind nebulae, SNR: Supernova remnants, unID: TeV uniden- 
tified sources. 

d : (1) Willingale et al. (2001); (2) Aharonian et al. (2004a); (3) DeLaney et 
al. (2006); (4) Aharonian et al. (2005b); (5) Aharonian et al. (2005c); (6) Ozaki 
& Koyama (1998); (7) Aharonian et al. (2005d); (8) Aharonian et al. (2004b); 
(9) Slane et al. (2001); (10) Brogan et al. (2005); (11) Aharonian et al. (2006); 
(12) Aharonian et al. (2005f); (13) Mukherjee & Halpern (2005); 
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Fig. 3. Spectra of the Src 1 (left) and Src 2 (right). Gray and Black crosses represent XISO+2+3 (added 
just for the clear viewing) and XIS1 data, respectively. The best-fit power-law models are shown with 
solid lines. Lower panels in the figures are data residuals from the best-fit models. 
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